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Abstract: Recently, an exact conformal mapping between soft gluons emitted from jets 
at large angle in e + e~-annihilation and those in the BFKL evolution of a high energy 
hadron has been proposed. We elucidate some remarkable aspects of this correspondence 
and use them to analytically compute the distribution and correlation of gluons in the 
interjet region. We also establish the timelike counterpart of Mueller's dipole model and 
discuss the resulting linear and nonlinear evolution equations. 
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1. Introduction 

Electron-positron (e + e~) annihilation into hadrons is one of the most well-studied high 
energy reactions that offers a broad arena for testing perturbative QCD predictions [1]. 
A number of fixed-order and resummed calculations with ever increasing precision have 
been developed for a variety of observables ranging from jet cross sections to event shape 
variables. The impressive agreement between these predictions and experiment witnessed 
over the past few decades undoubtedly represents a major success of perturbative QCD. 

While a large fraction of the theoretical activity in e + e~-annihilation has been centered 
around jet-related observables, there is a great deal of physics to be explored in regions 
between jets. A primary example is the energy flow E ou t [2-6], the total amount of energy 
radiated into a specified angular region away from the hard jets. The underlying partonic 
process that pertains to interjet observables is the multiple emission of soft gluons at wide 
angle. In perturbative calculations, large logarithms, called 'non-global logarithms' [7-15], 
of the type (a s lnQ/Ai) n appear, where Q is the center-of-mass energy and M. is a second 
hard scale Q > M > Aqcd characterizing the observable of interest. [M. = E out in 
the case of energy flow.] These logarithms arise due to emissions from secondary gluons 
and are therefore sensitive to complicated multi-gluon configurations in interjet regions. 
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Intertwined with the Sudakov logarithms [6] associated with emissions from the primary 
hard partons, their resummation is a challenging task in perturbation theory (see, however, 

[11])- 

Initially, the resummation of non-global logarithms to all orders was done numerically, 
in Monte Carlo simulations in the large- N c limit [7,8,15]. On the other hand, the authors 
of [9, 16] have succeeded in resumming logarithms in the form of evolution equations. 
Very interestingly, their results bear a striking resemblance to the BFKL [17, 18] and the 
BK [19,20] equations which have been hitherto discussed exclusively in the context of high 
energy (Regge) scattering. Indeed, the equations in [9] and [16] are almost identical in 
form to the BK and BFKL equations respectively, after merely replacing the kernel of the 
evolution equations as 



d 2 $l c 1 — cos 6 a b d 2 x c x 2 

4tt (1 - cos6> ac )(l - cos8 cb ) 2tt xl c x cb 



J ab 



2 



(1.1) 



Here, the left-hand-side is the well-known radiation function of a soft gluon (labeled c) 
from parent partons a and b (0 a b, etc. are relative angles between momenta.), whereas 
the right-hand-side is the dipole splitting probability in impact parameter space (a;'s are 
two-dimensional vectors and x a \, = x a — x^) which is the fundamental building block of 
the dipole formulation of the BFKL equation [21]. 

Such a resemblance naturally prompts one to seek a possible relationship between the 
two processes (e + e~ vs. Regge) at a fundamental level. In [22], it was pointed out that the 
map (1.1) is a conformal transformation, known as the stereographic projection. Moreover, 
in the strong coupling limit of M = 4 supersymmetric Yang-Mills (SYM) theory the same 
transformation exactly relates the final state in e + e~ -annihilation and the high energy 
hadronic wavefunction in impact parameter space. The way this latter result was derived 
(see [22] for the details) emphasizes that in the soft sector the two processes are one and 
the same phenomenon, the only difference being the choice of the coordinate system in 
which to express its physics content. In view of this, the correspondence (1.1) at weak 
coupling is hardly accidental, but must have a deep geometrical origin that goes beyond 
the perturbative framework. 

The purpose of this paper is to establish a detailed dictionary of the transformation 
rules and explore its physical consequences. The practical advantage in doing so is that on 
the BFKL side a number of exact analytical results are known with the help of conformal 
(or rather, the Mobius) symmetry. Making the most of the dictionary, one can obtain 
analytical insights which incorporate the effects of the resummation into the partonic final 
state in e + e _ -annihilation. We start by in the next section reviewing some basic facts 
about the stereographic projection which realizes the map ( |1.1| ). We clarify how the map 
correctly accounts for the subtle difference in kinematics between the two processes. Then 
in Section 3, we discuss the single and double gluon angular distributions and related 
observables in the interjet region. In Section 4, we construct the exact timelike analog 
of Mueller's dipole model using the generating functional technique. In light of this, the 
nonlinear evolution equations in the timelike and spacelike contexts can be treated in a 
unified fashion. Then in Section 5, we study the correlation of dipoles (heavy-quark pairs) 
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in the interjet region based on the results in [23,24]. Finally in Section 6, we look into the 
small-angle (distance) limit of the energy correlation functions in the dipole model and 
reproduce the OPE result in [25] . 



2. Stereographic projection 

The stereographic projection is a mapping between the unit sphere with coordinates 17 = 
(9,4>) and the two-dimensional plane x = (x 1 ^ 2 ). It is defined by the relations 

1 _ sin 9 cos (f) 2 sin sin <p 
X = l + cos# ' X = 1 + cosfl ' ( } 



or equivalently, 



2 1 I oc^ 1 

cos 9 = 7T, sin# = 7T, coscj)=- — r , sin0=- — -. (2-2) 

1 + x z l + x z \x\ \x\ 



The squared length transforms as 



( dx f = ITT ^2 ( d ° 2 + sin2 ed ^) = ^2 dn2 > ( 2 - 3 ) 

(1 + cos 9) z (1 + cos 0) z 



and the area element as 



d 2 Q = (I + cos 9) 2 d 2 x = ^d 2 x. (2.4) 

(l + x z y 



If one thinks of the sphere as being embedded in a three dimensional space (y 1 , y 2 , y 3 ), the 
stereographic map can be viewed as a part of the following conformal transformation in 
four-dimensions [26]; 

y -;tt> x = ' ( 2 - 5 ) 



2y+' y 2y+' ^2y 

where y = (y l ,y 2 ) and x ± = (x° ± x 3 )/\/2, y ± = (y° ± y 3 )/\/~2. 

We shall regard x as the transverse plane perpendicular to the direction of a high 
energy hadron. 1 The operator which measures the energy density at x is given by 



1 f°° 

£(x) = -j=j dx~T__(x+ = 0,x~,x) . 



(2.6) 



On the other hand, the coordinates Q = (9, 4>) are identified with the polar coordinates of 
quarks and gluons in the final state of e + e~-annihilation. The total four-momentum as 
measured in the y-coordinates is related to £ (x) via the following rules [25] 

P + = V2 J d 2 x£(x), (2.7) 
P~ = V2 J d 2 xx 2 £(x), (2.8) 
P = 2 j d 2 xx£{x) . (2.9) 



1 In the presence of conformal symmetry, it is natural to define it to be a dimensionless variable. When 
necessary, one can easily restore the length scale in the problem. 
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In particular the energy is given by 

E = J d 2 n£(n) = ^=(P + + P-) = j d 2 x(l + x 2 )£(x) , (2.10) 
where the operator [2,3] 

f°° 2 

S(n) = lim r 2 / dy° ni T *(y°,rn) = — m S{x) , (2.11) 

r->oo J Q (1 + COS 9)^ 

with ft = (sin 9 cos 4>, sin 9 sin <f>, cos 0) and r = \/y'( + y 2 , + y 2 ; measures the total energy 
flowing into the direction Q. In (|2.11| ), the factor of 2 accounts for the fact that to a high 
energy hadron with four-momentum (E, 0, 0, E) in the x-coordinates corresponds a virtual 
static photon with four-momentum [Q = 2E, 0, 0, 0) in the y-coordinates. 

In the following we shall be often interested in the energy and momentum of individual 
gluons rather than the total four-momentum. The transformation rules ( |2 . 7| ) — ( pT9| ) instruct 
us to make the identifications 

p^2xk°, (2.12) 
p° (1 + x 2 )k° , (2.13) 

where we employ the convention that the four-momentum of gluons in the timelike cas- 
cades (the y-coordinates) is denoted by while that in the spacelike cascades (the x- 
coordinates) is denoted by k^. [In fact, in the present approach the transverse momentum 
in the spacelike problem k does not appear explicitly, but only implicitly as the inverse 
of the transverse coordinates, x ~ k .] Note that equations ( 2,12j ) and ( [2.13 ) are not 
independent of each other due to the identity 

M = sin # = J^. (2.14) 
p u 1 + x l 



So far, the stereographic projection ( |2.1| ) has been introduced merely as a rule to 
associate particles living in two different coordinate systems. However, it turns out that 
this correspondence is preserved by the QCD evolution in the soft approximation. As 
observed in [22], the differential probability of emitting a soft gluon from a dipole (a 
quark-antiquark pair) with opening angle 9 a b 

_ d 2 n c l-cosfl afe _ 2 _ 

a s — ; — 7; 7— tt; ^— r = a s d il c K ab (il c ), a s = a s N c /Tr, (2.15) 

47T (1 - cos 9 ac ){l - cos 9 ch ) 

is exactly mapped via the stereographic projection onto the gluon emission probability 
from a dipole with transverse size x a ^ = x a — x^ in the spacelike process, 

a s X ° X 2 a \ = a s d 2 x c K ab (x c ) . (2.16) 

x ac x cb 

[Thus we use the same notation K in both cases but from the argument of K it should 
always be obvious whether we mean ( 2.15| ) or ( 2,16| ),] This implies that, at least to leading 
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logarithmic accuracy (and in fact in the large- iV c approximation), the high energy (small 
Bjorken-x) QCD evolution in the transverse plane is equivalent to the small Feynman-x 
structure of the interjet parton shower in e + e~ annihilation. In order to genuinely establish 
this statement and also to discuss its limitations we must, however, take a closer look into 
the details of kinematics. Indeed, there is a subtle but conceptually important difference 
between the two processes which needs to be addressed: When describing parton cascades 
in e + e~ annihilation, it is common to use the variable \p\ as the evolution parameter: 
One usually starts with a large value of \p\ set by the splitting of the photon into the 
quark-antiquark pair, and evolves the cascade towards smaller values of p with strong 
ordering in their magnitude \p\. It is then natural to use Yt = ln(E/\p\) as the evolution 
"time". In contrast, in the spacelike parton cascade of a high energy hadron the transverse 
momentum is more or less constant while one has strong ordering in energy, which in turn 
implies strong ordering in angle. One can then use either k° or the emission angle 9 as the 
evolution parameter with the evolution time Y s = ln(E/k°) or Y s = In 9, and to leading 
order the two choices should be equivalent. This difference in kinematics is an unavoidable 
feature of the multiple soft gluon emission in each case [16], and might be regarded as 
an obstacle against any attempt to find an exact mapping between the two processes. 
Remarkably, however, the stereographic projection automatically converts the nature of 
the evolution parameter into the desired form. Indeed, using ( |2.12 ) we have 



Due to our identification of x as the inverse transverse momentum k , the product 
\x\k° ~ k° /\k\ = 1/9 is indeed a measure of the emission angle, and thus we are lead 
to the identification Yt Y s . This assures that, to leading logarithmic accuracy, the cor- 
respondence works perfectly including the details of kinematics, and the whole machinery 
developed for the dipole formulation of the BFKL evolution can be used to analyze interjet 
observables in e + e~ annihilation, or vice versa. 

Let us conclude this section with a few additional remarks: 

It seems that the correspondence crucially relies on conformal symmetry, and as such, 
it may not hold, or at least needs to be modified, in the next -to-leading logarithmic (NLL) 
approximation in QCD where the running coupling effect breaks conformal symmetry. 
Indeed, as already pointed out in [16] the argument of the running coupling should be |p| 
and \k\ in the two cases, respectively. This does not agree with the rule ( |2.12| ) derived 
from a consideration of conformal symmetry alone. On the other hand, in J\f = 4 SYM 
theory which is conformal, there is a good possibility that the correspondence holds to all 
orders in the soft approximation, as indicated by the fact that it holds exactly in the strong 
coupling limit [22]. The recent NLL result reported in [27] is very encouraging from this 
point of view. In Appendix A we transcribe their result to obtain the timelike NLL dipole 
kernel. 

It is worth mentioning that the collinear singularity \p\ ~ 9 —* in the timelike emission 
kernel (glp maps onto the ultraviolet singularity |jb| ~ l/|fc| — > of the spacelike emission 



kernel ( 2.16[ ). Conversely, the ultraviolet region \p\ — * oo maps onto the collinear region 
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|fe| — * 0. Also note that in the timelike case this collinear singularity is responsible for 
generating angular-ordered gluons surrounding the primary quark and antiquark which 
eventually materialize into observed jets of hadrons with the multiplicity given by the 
standard double-logarithmic formula [28]. These collinear gluons are also mapped onto 
the transverse plane of a high energy hadron via the stereographic projection, although in 
this latter case they are irrelevant fluctuations which are basically invisible in the scattering 
process. In Appendix B, as another interesting aspect of this correspondence, we describe 
the interplay between the boost in the y^-frame and the scale transformation (dilatation) 
in the x^-frame originally noted in [25]. 



3. Interjet gluon distribution 
3.1 Single gluon distribution 

As a concrete example of the above mapping, let us compute the distribution of a single 
gluon emitted from a color dipole (a qq pair). We begin with the x-space and denote by 
£C a (6) the coordinates of the quark (antiquark). The single gluon distribution is 

d 2 N _ „ , , rdk° _ „ , , fHm) 

d 2 x 

where Y s = In 2 ^ k o and A is a cutoff. 2 Applying the stereographic projection, we obtain 



f dk° r\n(E/A) 

a s K ab (x) y — = a s K ab (x) J dY s , (3.1) 



d 2 N 1 d 2 N r HE/A) 



a s K ab (n) / dY t 
Jo 



d 2 n ~ (l+cos0) 2 d 2 x -"»"»v-vy o 

= *J-F& (32) 
4vr sin 2 9 J ' A 2 p 2 ' 1 ' 

where in the last equality we consider the back-to-back jets configuration 6 a = 0, b = ir, 
which corresponds to the choice x a = 0, x b = oo, see figure |l[ Note that we have made 
the substitution Y s -> Y t = Ya{E/\p\) following (|2.17|) . 

The energy distribution can be similarly computed. In the spacelike case we have 

£(x) = a s K ab (x) / k° dY s = a s K ab {x) / —— dY s . (3.3) 

Jo Jo z \ x \ 



Using ( 2.11 ), we find, for the back-to-back jets, 

2 

(1 + COS) 

a s (l + cos#) 2 f Ee~ Yt 



£(Sl) = Tr ——£{x) 



-dY t 



(1 + cos60 3 2vr sin 2 6» J 2 sin 0/(1 + cos 
™* 1 f d P 2l -4, (3.4) 



47r sin 3 9 J p 2 

in agreement with [5]. We note, incidentally, that in the x space the energy distribution 
far away from the parent dipole falls as ~ l/|a?| 5 as can be seen in (|3.3j). 



2 The limits on the Y s integral would imply 1 > 6 > A/E in terms of the emission angle 8 w |fe|/fe° 
A/|x|fc°. 



-6- 



Figure 1: Stereographic map between the unit sphere and the transverse plane. 



3.2 Resummation to all orders 



The results (3.1) and fl3.2|) can be viewed as the leading order term in the expansion in 
powers of a s Y. 3 If a s Y becomes of order unity, one has to resum all the higher order 
terms (a s Y) n consistently. In the spacelike case, this can be done in Mueller's dipole 
model [21,29]. In this approach the average total number of dipoles (xi,X2) within a 
rapidity interval Y contained in the parent dipole (x a ,x&) is given by 4 

ny(x ab ,x 12 ) = — \- V * U + 4) e &sXMY 



X 



2^ 2 ) 2 V^oo(2^H 4 
J cPxeE 1 -^^, x 2c )E h ~ h (x ac , x bc ) , (3.5) 



where x(n, i/) is the usual BFKL eigenvalue, and E ' is the eigenfunction of the SL(2,C) 
Casimir operator (see, e.g., [30]) with h = + iu, h = + iv. 

Formally, the total number of dipoles (equal to the total number of gluons plus one) 
is given by the integration 

N Y = j d 2 xid 2 X2n Y {x ab ,xi2) . (3.6) 

As we shall soon see, the integral is actually divergent. Ignoring this fact for the moment, let 
us analyze its structure from the viewpoint of conformal symmetry. The d 2 x c integration 
in ( |3.5| ) gives a function (see below) only of the anharmonic ratio 

ZabZl2 ( taI1 ~ taI1 ^ e ^) ( taI1 l"^ 01 ~ taI1 9 -2 ei4 / 

P ~~ = ~7 \ / 7~ j (3.7) 

z a \Zb2 /tan - tan J ( tan fe^ - tan ^e^ 2 



3 In the following we do not distinguish Y t from Y s . 

4 Our normalization of ny differs from the usual one in the literature by a factor of 1 / (2-KXi 2 ) ■ Note also 
that despite our somewhat sloppy notation, ny in fact depends separately on x a , Xb, &i and The same 
remark applies to the other distributions to be defined later. 
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and its complex conjugate. [We have here introduced the complex coordinates z a = x\+ix 2 a , 
etc.] The square of this ratio is 

I i2 _ x lb x i2 (1 -cosg afe )(l -cosgi 2 ) . 
lPl x 2 al x 2 2 (l-cos0 al )(l-cosM' K ' } 

Moreover, the integration measure 

d 2 xid 2 x 2 _ d 2 nid 2 tt 2 . . 

(x 2 2 ) 2 ~ 4(l-cos0 12 ) 2 ' 

is conformally invariant so that Ny can be expressed in a manifestly conformally invariant 
way. Keeping only the n = term which is dominant at high energy, we get 

d 2 Xid 2 x 2 n Y (x ab ,x 12 ) = / . 2 . „ f(p, p) , 

J { X 12) 

= [ d 2 n,d 2 n 2 n Y (n ab , n u ) = \ [ ,f nid J\ 2 f( P , P) , (3.10) 

J 4 J (1- cos e 12 y 

where [30,31] 

/oo . 
(fa/ e a ^W y (Mp| 2 ~ 27 2 iq(l - 7, 1 - 7 , 2 - 2 7 ; p) 2 Fi(l - 7, 1 - 7, 2 - 2 7 , p) 
-00 

+bt\p\ 2 ^ 2 Fi( 7 , 7,27^)2^(7, 7 ,2 7 ;p)) , (3.11) 

with 7 = i + ii/ = /i, and 

= ^2 4 ^ £(| - it/)r(l + iv) 

K = T(\+iu)T(l-iu)- (3 ' 12) 

In the second line of ( |3.1C| ) we are naturally led to define the timelike analog of the single 
dipole distribution ny (£l a b, ^12), that is, the total number of dipoles with opening angle 
O12 contained in the parent dipole tt ab within a rapidity interval Y . A related distribution 
(integrated over (Q,\ + £l 2 )/2) was previously introduced in [16] in the context of the heavy 
quark pair production in e + e~-annihilation. We shall discuss more about this in the next 
section. 

The generalization of the angular distribution (|3.2j ) would be 

d 2 Ny (1 + X 2 ) 2 d 2 Ny (1 + X 2 ) 2 



2 J d 2 x 2 ny(x ab ,xi 2 ) . (3.13) 



(PQi 4 d 2 x 1 4 

If we were to assume the d 2 x 2 (or d 2 £l 2 ) integral to be convergent, then from dimensional 
analysis the result would be proportional to 

/ d 2 x 2 n Y (x ab , x 12 ) oc —» , (3.14) 
J xf 

in the limits x a — * 0, x b —* 00, and therefore, 

^oc<I±# 1 . (3.15) 
d 2 Qi Ax{ sm 2 6»i 
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Figure 2: The angular correlation of gluons in a highly boosted frame of the qq pair (left) and in 
the center-of-mass frame (right). 



This has the same angular dependence l/sin 2 # as in the lowest order result ( |3.2| ), and in 
fact it is the unique possibility consistent with boost invariance [5]. However, as already 
mentioned the integral is divergent due to the singularity at x 2 = X\ (or Q% = Indeed, 
when ~ Q\ 2 is very small, \p\ is much smaller than unity. We may then approximate 
as 2^ ? i(---! p) ~ 1 and find 

!Nh~J-J dn \i- C ose 12 )^ - (3 - 16) 

For small values of #12, the anomalous dimension 7 at the saddle point is between and |, 
hence the singularity is not integrable. It would be interesting to see whether this problem 
is cured in a more refined treatment including higher order corrections (in particular, the 
energy conservation) to the BFKL approximation. 

On the other hand, the doubly-differential gluon distribution is finite and can be 
identified with the dipole density itself. 5 When \p\ <C 1, a simple analytical estimate is 
possible. This encompasses two physical situations (see fig. ||) which are mathematically 
equivalent: (i) Radiation from a highly boosted qq pair. In this case 8 a b <C 1 and 

| |2 W ^(1-C0SM 

lPl ~2(l-cosM(l-cosM { V 

Evaluating the 7 integral in the saddle point approximation, we find 

d 4 N Y e SsX(7s)Y 



d?VLid 2 VL 2 (l-cos6'i2) 1+ T s (l-cos0 a i) 1 -^(l-cos0 a 2) 1 -^ ' 



(3.18) 



Some care must be taken in this identification since the two gluons are not selected randomly but are 
constrained such that they are neighboring in the color space. Though we presume that such a concern is 
immaterial in the large N c limit at the level of the two-gluon distribution, the following results may admit 
somewhat different interpretations such as the measure of color flow at large angle. [Note that even the 
softest gluons carry color.] 



-9- 



where the saddle point < j s < \ is determined by the equation 

As an example of the value of j s , we note that 7 S = 0.47 for \p\ 2 = 0.1 and a s Y = 2. As 
\p\ 2 decreases further the saddle point moves slowly towards lower values. 

(ii) The back-to-back jets case with 9 12 <C 1. Setting 9 a = 0, 0& = ir and therefore, 

< 1 , (3.20) 



sin c/i 
we obtain 



(3.21) 



where the saddle point is again given by (3.19). 



When 9\2 becomes large, \p\ reaches values of order unity. In this regime the hyper- 
geometric function 2F\{...; p) has to be fully retained. As an example, let us consider the 
azimuthal correlation of two gluons around the jet axis by taking $l a = (0,0), = (vr,0), 
ill = (f ,0) and 0,2 = (f ,<!>)■ Equation ( ]3 . 7| ) becomes 



j0 = l_ e # ) |p| 2 = 2(1 - cos0) . (3.22) 

This leads to 

^^ 2 =Aj^^^2F l {l- 1 ,l-^2-2r, P )2F l (l- 1 A- 1 ,2-2 1 ,p). 

(3.23) 

We have numerically integrated the right-hand-side of fl3.23|) . The result is shown in fig. || 
as a function of <f> for two different values of a s Y. We see that the strong correlation in 
the collinear direction 9\2 = <f> <C 1 as described by ( |3.21| ) decreases towards the backward 
direction and eventually reaches a minimum at <fi = ir. 

Another interesting physical quantity is the (pseudo-)rapidity correlator of gluons in 
the interjet region. The rapidity here is defined as 

rj = In cot-, (3.24) 
in terms of which the anharmonic ratio reads 

p = 1 - cot y tan ?l e i U»-<M = 1 - e V* . (3.25) 
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Figure 3: The azimuthal correlation of two gluons around the jet axis for ct s Y = 2 (upper curve) 
and a s Y — 1.5 (lower curve). The horizontal axis is <f> = 9\2- 



The appearance of the relative rapidity rj = r/i — r]2 is a consequence of boost invariance [5] . 
Integrating ( |3.23| ) over the azimuthal angle, we get 6 



d*N- 



y 



8vr 



2- 



-CO \P\^ 



d cos 6*i d cos 82 jo 

x 2 F 1 (1 - 7 , 1 - 7 , 2 - 2 7 ; p) 2 F 1 (l - 7, 1 - 7, 2 - 2 7 , p) . (3.27) 

The large 77 behavior of the above integral can be estimated as follows. We use the identity 

2 2a 



2 Fi(a, a, 2a, p) 



T(a + \/2){-p)- a Ql_ 1 (l - 2/p), p i (0, 1) 



\/vFr(a) 



(3.28) 



where Qu is the Legendre function of the second kind. Since \p\ « e 11 3> 1, we can neglect 
the p dependence of the Legendre function. [Strictly speaking this function goes to an 
infinite constant in the p — > 00 limit.] Using a = 1 — 7 in our case, we get 



d 4 N Y 



dcos Videos 9 2 



d^y 



1 



| p |2+2 7 | p |2-2 7 



1 



4 x I(Y) 



(3.29) 



where 1(1") is the value of the 77-independent integral. We thus see that at large rj the 
correlator decays as e _4r? . Fig. ^| is the result of a numerical integration of ( 3.27 ) as a 
function of r]. The correlator indeed shows an exponential decay e _cr? with c ~ 4 already 
when rj > 1, more or less independently of the value of a s Y. In comparison, we note that 
the energy correlation function in the lowest order (two-gluon) approximation exhibits a 
similar exponential decay with c = 3 [5]. 



6 It is physically obvious that the correlator should be invariant under the sign flip 77 — > —r\, thus it is a 
function only of \r\\, or rather, e v + e~ v = 2cosh?7. In (3.27), this can be seen by using an identity of the 
hypergeometric function 



2 F 1 (a, A 7, «) = (!- z)- a 2 Fx a, 7-/3, 7 



(3.26) 



or more easily by choosing Q a = ("",0), Qb = (0,0). 
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Figure 4: Logarithm of the right-hand-side of ( 3.2? ) for a s Y = 2 as a function of r] 



4. Equivalence of the dipole formulations 

The single dipole distribution ny(O a b, Q12) denned and studied in the previous section is 
the exact timelike counterpart of the spacelike distribution ny (a; a &, X12). In this section we 
show that this correspondence can be generalized to the complete equivalence of the dipole 
formulation of the timelike and spacelike cascades. The generating functional in Mueller's 
dipole model is given by 

00 „ fn— 1 \ 

Z ab [u,Y} = ^2 \Y[d 2 Xi \P n (x 1 ,- --Xn-i ;Y)u 1 u 2 --- u n , (4.1) 

n=l^ \i=l / 

where Ui = u(xi-i,Xi) is an arbitrary ' weight' function for the i— "th dipole and {xq^Xji) = 
(x a ,Xb) is the parent dipole. The function P n is the probability distribution to have n 
dipoles in a cascade evolved up to rapidity Y. From probability conservation, Z[u = 1] = 1. 
The equation satisfied by P n is 

n „ 

d Y Pn(xi,---x n - 1 ;Y) = -a s ^2 / d 2 z K i ^ i (z)P n (x 1 , ■ ■ ■ , x n -r, Y) 

i=l •* 

n-l 

+a s ^2 Ki-x^+iix^Pn-xix!, ■ ■ ■ Xi-i, x i+1 , ■ ■ ■ x n -x;Y) . (4.2) 



This has a simple physical interpretation as a gain-loss type of equation. The first term 
on the right hand side is the "loss" term and describes the total probability for the dipole 
configuration to disappear via all available decay channels. The second term is the "gain" 
term which is the probability to obtain the given n-dipole configuration from all possible 
n — 1 dipole configurations. The single dipole density is calculated from Z a b via 



5Z ab [u,Y] 

n Y {x ab ,x 12 ) - 



Su(xi,X2) 



(4.3) 

u=l 
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and it obeys the BFKL equation 

d Y n Y {x ab ,xi 2 ) = a s J d 2 x c K ab (x c ) [n Y (x ac , x 12 ) + n Y (x cb , x u ) - n Y (x ab , x 12 )] 

= a s J d 2 x c K ab (x c ) ®n Y (x ab ,x 12 ) . (4.4) 
On the other hand, the evolution equation for the generating functional is 7 

d Y Z ab = a s j d 2 x c K ab (x c )(-Z ab + Z ac Z cb ) . (4.5) 



In particular, if one chooses u(xi_i, Xi) = s(xi^i, Xi), the dipole S'-matrix, ( f4.5| ) is nothing 
but the Balitsky-Kovchegov (BK) equation [19,20] for the total S'-matrix Z ab = S ab . The 
BFKL limit of the BK equation is obtained by defining the T-matrix t = 1 — s, T = 1 — S, 
and expanding T to linear order in t. One finds 

d Y T Y (xab) = a s J d 2 x c K ab {x c ) [T Y {x ac ) + T Y (x cb ) - T Y (x ab )] , (4.6) 

with 

T Y (x ab ) = J d 2 x 1 d 2 x 2 t(x 1 ,x 2 )n Y (x ab ,x 12 ) . (4.7) 

Let us now turn to the timelike cascade and construct the probabilistic interpretation 
of the dipole evolution. First we recall the angular part of the real emission probability of 
k gluons in the soft kinematics [28], 

dV k (p a ,Pi,---Pk,Pb) = d 2 n l ---d 2 n k - 1 — C °! 6ab — (4.8) 

(1 - COS# a l)(l - COS#l 2 ) •••(!- COS& kb ) 

In the dipole language this can be viewed as the production probability of (k + l)-dipoles 
(a, 1), (1, 2), . . . , (k, b) obtained by multiplying together successive factors of the kernel Kij 
according to a particular history of the dipole cascade. In the final result, the dependence 
on all the intermediately formed dipoles disappear and each factor in the denominator 
associates with one final dipole, while the numerator comes from the decay of the original 
dipole (a, b). Under the stereographic projection, (|4.8j) becomes 



x 2 

dV k (x a , xi, ■■ ■ , x k x b ) = d 2 xi ■ ■ ■ d 2 x k — — ^ T , (4.9) 

^01^12 ' ' ' X kb 

which is precisely the corresponding result for the spacelike cascade. Here the denominator 
contains a product of the squared length of each final dipole in the cascade, and the 
numerator contains the squared length of the original dipole. 

One can similarly show that the virtual contributions are also mapped onto each other 
via the stereographic projection since they are simply obtained by integrating the real 
emission kernel K over final state coordinates. Therefore one is guaranteed to have the 



See [32] for an illuminating discussion on the consistency between (4.2) and ( 



- 13 - 



exact timelike analog of the generating functional techniques described above. This is 
readily achieved by applying the stereographic projection to ( |4,1| ) 



oo ~ /n— 1 J2q. \ 
Z ab [u,Y] = Y,l II (i + cos V)2 ^(^i' • • • *>i"2 • • • n " • ( 41 °) 

n=l \i=l / 



It is natural to define the probability distributions in the timelike case as 

nr^d+ooe 



With this definition one can easily check that the evolution equations satisfied by P n and 
Z ab in the timelike case are identical in form to the corresponding equations (fl^) and ( [4.5| ) 
after replacing a?j — ► everywhere. 

For certain applications, mainly in the timelike context, the need arises to specify the 
energy (or rapidity) of individual partons. Equation (4.1) is not suitable for such purposes 
since P n (...,Y) is already integrated over energy. The more appropriate definition of the 
generating functional would be 8 

°° n ( f f E dk°\ 
Z ab [ Ul E} = J2Y[{ d 2 Xi -j±)P n (x 1 ,---,x n )u 1 u 2 ---u n , (4.12) 

n=l i=l ^ 1 ' 

where we have introduced the probability distribution of n gluons (rather than dipoles) 
which are ordered in energy, and accordingly, let the source function depend on one coor- 
dinate as well as on energy 

m(xi-x,Xi) u(xi,J$) . (4.13) 

The equation for Z is [21] 

Ed E Z ab = a s j d 2 x c K ab (x c )(-Z ab + u(x c , E) Z ac Z cb ) . (4.14) 

In the timelike case, the energy integral becomes dk° /k° = dY s — ► dYt = dp°/p° and, 
due to the correspondence ([Q| )- (|4~9|) , it follows that \\ i dx i P n ({x}) — > Y[ i dQiP n ({Q}). 
Therefore, the equation for Z in the timelike case is identical to (|4.14[) except that x c is 
replaced by f2 c . 

The energy flow observable considered in [9] is the probability that the total amount 
of energy emitted into a specified interjet region C ou t is less than E ou t <C E. The nonlinear 
evolution equation derived there has precisely the structure ( |4.14| ) with the weight function 

(SU,Pi) = ©m(a) + e- p °/^G oui (a) . (4.15) 



u 

This follows from the kinematical constraint 

n-1 



@(E out - Pfi ~ II ( -(^) + e~ P * /Eout ®out(ni)) , (4.16) 

j&Cout 1=1 



3 In fact, this was the original definition employed in [21]. 
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where Q ut(ty is the support function nonvanishing in C ou t and Oj n (f2) = 1 — @ out (£l). 

More generally, for each physically motivated choice of u(Qi,p%) or -u(f2j_i, fy) (u(xi, A:f ) 
or u(xi-i,Xi) in the spacelike case) the generating functional becomes an observable. Al- 
though the interpretations of the resulting evolution equations (with or without a u-factor 
in the nonlinear term) vary drastically depending on the context, mathematically they are 
equivalent and can be mapped to one another via the stereographic projection. 

5. Dipole pair density 

Having established the timelike version of Mueller's dipole model, we can now in principle 
study arbitrary higher order correlations among dipoles in the interjet region of e + e~ 
annihilation. Physically this is relevant to the number correlation of heavy-quark pairs. 
Similarly to the single dipole distribution ny{x ab ^x c i) fl4.3|), one can define the fc-dipole 
inclusive distribution by differentiating k times the generating functional 



(*)/ v 1 s k z ab 

n Y {x a b;xiir,X22', -,Xkk') - — 



• (5-1) 

u=l 



k\ 5u(xi, Xi/)5u(x2, X2>)...5u(xk,Xk') 

The corresponding distribution in the timelike case can be immediately inferred from ( 4.1l| ). 
In the case of k = 2 ('the dipole pair density' [29]), we find 

(2) 1 (2) 

n\r>(Q ab ,Q lv ,n22>) = TT7 WT2 U Y ( x ab,X n >,X22>) • (5.2) 

11^11'22'V 1 + cos Vi) 

(2) 

In [33], an exact integral representation of n Y (#ab> ^ii'j ^22') has been derived. [See also 
[34].] Unfortunately, the expression is too complicated to be evaluated in full generality. 
However, in certain limits analytical results are available [23,24]. These include the large 
parent limit 

laJobl > \x\2\ > \x\v \1\x22' I , (5.3) 

and the small parent limit 

\x a l\, \x a2 \, \x 12 \ > |£C a6 |, |a?ll'|j \ X 22'\ ■ (5.4) 

As observed in [23], the above two configurations are transformed to each other via a 
conformal transformation. This means that the results can be unified in a single expression 
which involves anharmonic ratios 

J/ ll' x 22' \ x 12 / \ J 'al Jj b2/ 

where the anomalous dimensions 0<7 s <7^<1/2 are determined from certain saddle 
point conditions. 9 Using (|5.2[), ([5.5[) and the stereographic projection, we find 



(2), n o v e 2xW)y /(l-cos^nOa-cos^U 1 -^ 

ny{n ab ,i} 11/ ,n 2 2')- ' 



(l-COS0 n /) 2 (l-COS022') 2 V (l-COS^l 2 ) 2 



1-7* 



' (!_- cos6> ab )(l - cos 6*i 2 ) . . 
l i (l-cos0 al )(l-cos^ 2 ), " J ' 



9 See [23] for details. Our normalization of n' 2 - 1 differs from that in [23] by a factor (pu\ v x\ 2 ^ 1 . Also 
the anomalous dimension is redefined as 7 — > 1 — 7. 



- 15 - 



The case ( p^q ) includes the back-to-back configuration a b = tt. [See fig. g, but this time 
the gluons 1 and 2 are replaced by the dipoles 11' and 22'.] The expression in ( |5.6| ) then 
reduces to 



(2 ) e 2 *W) Y 1 



' (5-7) 



(l-cos^-^l-cos^) 1 - 7 

where the result is valid when 2(1 — cos #12) <C (1 — cos o i)(l — cos O^). On the other hand, 
the case ([TJ) corresponds to radiation from a highly boosted qq jets such that 9 a b -C 1. 
Equation (p\q ) then reduces to 

(2) e 2x(~f' s )Y 1 

n y (O a6 ,O 1V ,fi 22 ~ 2+ 27 ^2+2^ (1 _ cos012) l-2^ +7s 

/)2-2 7s 

afe (5.8) 



(l-cos^Oi-T^l-cos^a) 1 -^ ' 
It is interesting to compare ( |5.7| ) and ( |5.8| ) with the two-gluon correlation function ( 3,18j ), 



( 3.21 ). Putting aside a possible numerical difference in 7 S , we see that the growth of the 



correlation as 612 —* in (|5.8[) is significantly weaker than in ( 3.18j ). This is of course due 



to the color screening effect of dipoles. Moreover, when Qyi becomes comparable to either 
d\\> or #22' the correlation function converges to a finite value. This conclusion cannot 



be reached from (5.6) which assumes #12 3> #11', #22'- Rather, it follows from a proper 
(2) — 

evaluation of n Y in the regime #12 ~ #11', #22' (or | as 12 1 ~ l a; 22'|) as was done in [24]. 

[In the spacelike case this regime deserves special attention in association with the BK 
equation.] 

6. Energy correlation function 

In [25], the small angle limit of the energy-energy correlation functions in e + e~ annihilation 
was studied using the method of the operator product expansion (OPE) 

(g(ni)g(n 2 )) ~ ,^,^(3) , (012-O), (6.1) 

where is the anomalous dimension of the twist-two operators with spin j. Essentially 
the same OPE applies to the high energy hadron problem [22] 

(S( Xl )£(x 2 )) ~ ]xi J +2j{3) , (*i2-0). (6.2) 

It is straightforward to generalize these results to higher point correlation functions 

mm) • • • *w> - ^,-2+2,(^1) . ( 6 - 3 ) 
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and 



(£(x 1 )£(x 2 ) • ■■£{x k )) 



1 



(6.4) 



x 



2fc-2+27(fe+l) ' 



where we let the relative angles Oij and coordinates \xi — Xj\ all go to zero while keeping 
their ratios fixed, and denoted their representative values by 9 and x, respectively. 

The above results are valid not only in N = 4 SYM where they were originally derived 
[25], but also in any conformal theory. In QCD, they are expected to hold in approximations 
where the running of the coupling is neglected (as in the leading order BFKL). In this 
section we study the energy correlation functions in the dipole model and provide a concrete 
physical picture of how these singular behaviors are dynamically generated in the parton 
evolution. While there is the usual caveat about discussing energy-related observables in 
energy-nonconserving approximations, we believe that the dominant configuration which 
we shall identify below continues to be the relevant one even after a proper implementation 
of energy conservation both in the timelikc [35] and spacclikc [36] gluon cascades. 

6.1 The energy two point function 

We first consider the energy two-point function. Though the primary interest in the cor- 
relation functions arises in the context of e + e~ annihilation where they are experimentally 
measurable, we find that the actual calculations are a little more transparent in the x- 
space. 10 Of course, via the stereographic projection all the equations to follow have a 
timelike counterpart. We employ the energy ordering scheme and write Y = \n(E/k°). 
The energy two-point function can be generically written as 



where n(x\,Yi,x 2 , Y 2 ) is the number density of gluon pairs with rapidities Y\ and Y 2 located 
at x\ and x 2 , respectively. At high energy, we would like to relate this quantity to the 
single dipole density ny (x a b, x\ 2 ). [Henceforth we suppress the dependence on the parent 
dipole coordinates x a b in ny and set |x a j,| = 1.] 

In general, there is no universal relation between n(x\, Yi, x 2 , Y 2 ) and ny{x\ 2 ); the 
former carries information on the rapidity of the gluons 1 and 2 (in a cascaded evolved 
up to Y — > oo), whereas the latter simply counts the total number of dipoles {x\,x 2 ) in 
a cascade evolved up to Y without specifying the energy of their constituents. However, 
as long as the singular behavior in the limit x\ — > x 2 is concerned, the two quantities 
can be linked by the following argument: For a given gluon pair with rapidities Y\ and 
Y 2 such that Y 2 > Yi, the strongest correlation comes from the moment in the history of 
the cascade when the two gluons in question were realized as a dipole as a result of the 
emission of the gluon with rapidity Y 2 . The number of such pairs at the time of creation 

The schematic argument given in [22] overlooks the subtleties to be raised in the following. Our refined 
argument is physically more correct and readily generalizes to higher order correlation functions. 




(6.5) 
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is related to 9y 2 ny 2 (^12) where the derivative is because the gluon 2 has rapidity exactly 
equal to Y 2 (rather than integrated over rapidity). This is, however, still inclusive in the 
rapidity of the gluon 1. Therefore, we are led to define the "unintegrated" dipole density, 
n u {x\ , Y\ , x%, Y 2 ) , as follows 

dY 2 n Y2 (x 12 ) = / dY 1 n u (xi,Y 1 ,x 2 ,Y 2 ) . (6.6) 
J 

To proceed, we observe that n u is approximately independent of Y 2 . In order to see 
this, consider a simple toy model of the dipole cascade in which we neglect the transverse 
dimensions and assume each dipole to split after a fixed rapidity interval 11 Ay. At Y = 0, 
there is only one dipole composed of a quark and an antiquark both having Y = which 
we denote as (0,0). This dipole emits a gluon at rapidity Ay and splits into two dipoles 
(0, Ay) and (Ay, 0). In the second step these two child dipoles split into four dipoles 
(0,2Ay), (2Ay,Ay), (Ay,2Ay) and (2Ay,0). The process stops after n = Y tot /AY 
steps, with Y to t being the total rapidity, and generates 2 n = 2 Ytot ^ AY = exp{ (In 2/ Ay )Yt t} 
dipoles. One can easily see that there are 2^/ AY ^~ l gluons having rapidity Y. In particular, 
half of the total number of gluons have the maximal value Y = Ytot- This means that all the 
dipoles have the rapidity composition either {Ytot-, Y) or (Y, Ytot) and their number depends 
only on Y, and not on Ytot- I n the full model the rapidity distribution is more smeared 
out than in the toy model, but it will again be true that the density n u (xi, Yi, x 2 , Y 2 ) is 
largely determined by the gluons at Y±, since their number is exponentially smaller than 
the number of gluons at Y 2 in the presence of the strong energy ordering Y2 3> Y\. 

We thus write n u = n u (xi,x 2 , Y\) and get 

d^n Yl {xi2)=n u {x 1 ,x 2 ,Y 1 ), (6.8) 

as the number density of dipoles (xi,x 2 ) having the rapidity composition (y,y2) in a 
cascade evolved up to Y 2 . Later in the evolution (i.e., at Y > Y 2 ) these dipoles may or 
may not split, but the gluons 1 and 2 remain fixed at their original positions (xi,x 2 ). 
Since the total probability of all the decay channels is unity, the most singular contribution 
(as x\ 2 — ► 0) to the two-gluon distribution n(xi,Yi,x 2 ,Y 2 ) in (6.5) which is in principle 



defined in a cascade evolved up to Y = 00 is actually frozen at Y = Y 2 (or at whatever 
value Y > Y\) and is given by ( |6.8[ ). 

We now recall the integral expression ( |3.10D for ny(xi 2 ) in the limit x\ 2 — > 

nAXl2) = xjJ ^^0)' (6 ' 9) 
where the contour j-integral goes around the branch cut in the Jacobian c(j) 

dj ~ 2^140^(3)0' - jo 



c(j) oc^^ , = J=. =. (j = l + 4a s ln2) (6.10) 



11 In the full model one has 



where p is a small cutoff. 
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Using this, we find the energy-energy correlator 

rY 2 



poo pY 2 

(£{ Xl )£ (x 2 )> ~ 2E 2 / dY 2 e- Y2 dYte'^ d Yl n Yl (x 12 ) 
Jo Jo 

= 2 ^y ^ 3-, i, 12 i 2 ^)- (6 - n) 

Deforming the contour and picking up the pole at j = 3, we find 

(f (Xi)£ (x 2 )) ~ + -L / dj • • • ^-y . (6.12) 

Since < 7(3) for j > 3, the contribution from the second term is subdominant in the 
limit x% 2 — > 0. We therefore recover the OPE result (6.2). 

6.2 The energy three point function and beyond 

The energy three-point function is given by 

poo pYs pY2 

(£(xi)£(x 2 )£(x 3 ))=3\E 3 dY 3 e- Y * dY 2 e~ Y * dYie- Yl n( Xl ,Y 1: x 2 ,Y 2 ,x 3 ,Y 3 ) , 

Jo Jo Jo 

(6.13) 

where we n(...,x 3 ,Y 3 ) is the distribution function of three gluons. The leading singular 
behavior as X\j — * arises from the configuration where the three gluons are realized as two 
contiguous dipoles which, for Y 3 > Y 2 > Y\, are {x\,x 3 ) and (x 3 ,x 2 ). This occurs when the 
gluon 3 is emitted from dipoles (x\,x 2 ) whose number density is given by (|Bij|). Thus to 
get the number density of the dipole pairs (331,033) and (x 3 ,x 2 ) with rapidity composition 
(Yi, I3) and (Y 3 ,Y 2 ), we just need to multiply (|6.8|) by the emission probability. This gives 
the unintegrated density 

n u (x 1 ,Y 1 ,x 2 ,Y 2 ,x 3 ,Y 3 ) « a s K 12 (x 3 )d Yl n Yl (xi 2 ) . (6.14) 

We now substitute this quantity in (|6.13|) as a dipole model analog of the three gluon 
distribution. We remind the reader that this is justified only for the singular contribution 
in the limit Xy — ► 0. Writing Xij ~ x and therefore ^12(333) ~ -4-, ( |6.13| ) then becomes 

poo pY 3 pY 2 

(£( Xl )S(x 2 )S{x 3 )) ~ i? 3 / dy 3 e" y3 / dY 2 e~ Y * / (fYie^ 1 

jo Jo Jo 

a s E 3 f dj , (j - l) 2 1 



2m v '2(4 - j) |x| 2 ^0') 
"^otW (6 - 15) 



in agreement with (|6.4|). 
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It is now straightforward to extend the above result to the energy /c-point correlation 
function 



/•OO l>Y2 

{£{x{)£{x 2 )---£(x k )) = k\E k dY k e- Y *--- / dY 1 e- Yi n(x 1 ,Y 1 ,x 2 ,Y 2 ,...,x k ,Y k ) . 

Jo Jo 

(6.16) 

The leading singular contribution to the fe-gluon distribution n(..., x k , Y k ) with Y k > 
Y k -i > ■ ■ ■ > Y\ comes from the process where we start with the lowest rapidity pair 
{x\,x%) and successively emit the remaining gluons, creating first the pair (031,033) and 
(033,032), then emitting the gluon 4 from either of these two dipoles, and so on. One then 
gets the fully unintegrated distribution 12 

a k s ~ 2 K 12 (x 3 )K 3 2(x^---K k _ lt 2(x k )d Yl nY 1 (x 1 2) ~ (^§) ^ny^x^). (6.17) 



Substituting this in ( 6.16j ) as a proxy for the /c-gluon distribution n(...,x k ,Y k ), one finds 
the result 

1 roo rY2 

(£(xi)£(x 2 ) ■ ■ ■ £(x k )) ~ a k ~ 2 E k -^ k -^ ) J dY k e~ Yk ■ j dY x e~ Y ^ 

~ E k a k ~ 2 (6 18) 

s | £C |2 7 (fc+l)+2fc-2 • yv.10) 

Summarizing, in the dipole picture the most singular contribution to the fc-point en- 
ergy correlation function comes from the moment in the history of the dipole cascade when 
the k gluons are realized as a single chain of k — 1 dipoles. This chain represents succes- 
sive emissions of energy ordered gluons at comparable separations (collinear emissions at 
comparable angles in the timelike case). Later in the evolution the chain decays, but the 
k gluons remain in their original positions. Integrating over all possible values of rapidity 
at which this chain is formed, one recovers the OPE result. 

In strongly coupled J\f = 4 SYM the AdS / CFT correspondence allows one to represent 
the energy correlation function as the scattering amplitude between the photon (the pri- 
mary qq pair) and the graviton (energy operator insertions, or calorimeters) [25]. Its small 
angle limit is associated with the Regge behavior of the string S-matrix and is dominated 
by the t-channel exchange of a massive string state ('Pomeron') with an anomalous di- 
mension 7(1 + k). A similar interpretation is possible in QCD by identifying the £-channel 
object with the BFKL Pomeron. Naively, one may expect from the relevant configura- 
tions that the /c-point function would be related to the (A; — l)-ple dipole density riy ^ 
(cf. ( |5.1[ )) which involves the exchange of (k — 1) Pomerons with an anomalous dimension 

7 = x~ 1 (^(k ~ 1) m 2) [37]. However, specifying the rapidity of k gluons amounts to taking 

(ft— 1) 

the £;-th derivative of n-y with respect to Y. Each Y-derivative lowers the degree of 

12 This expression is valid for a particular color ordering of the dipoles, but since in the end we let all 
\xi — Xj\ go to zero with fixed ratios, the result holds for any ordering. 
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multiplicity, thus one has a cascade of distributions — ► n^ k ~ 2 ^ — > n^ k ~ 3 ^ • [38]. 
The fully unintegrated distribution obtained in this way involves the single dipole density 
or the single Pomeron exchange with an unusual value of the anomalous dimension 7(1 + k). 
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A. Next-to-leading order dipole kernel in e + e —annihilation 

Recently, Balitsky and Chirilli have derived the spacelike next-to-leading logarithmic 
(NLL) dipole kernel in M = 4 SYM which turned out to be conformally (Mobius) symmet- 
ric [27] (see, also, [39]). As we noted in Section 2, in this theory there is a rather strong 
indication that the stereographic projection works to all orders of the soft approximation. 
Assuming this to be correct, one can immediately obtain from their result the NLL kernel 
for the dipole evolution in the timelike case by applying the stereographic projection. The 
evolution equation for the dipole density which generalizes the leading order equation ( \i-4 ) 
is 



dyny(Q 



ah) 



a* 1 — a. 



71 

12 



d 2 Q. c K ab (tt c ) [n Y (VL ac ) + n Y (Vt cb ) - n Y (fl a b)] 



+a 2 s J d 2 n c d 2 n d K' ab (Q c ,n d )n Y (n cd ) , 



(A.l) 



where 



K b (n c ,n d ) = —5 



(i 



cos 



8tt 2 1 (1 - cos6> ac )(l - cos6> cd )(l - COS 6 db ) 

(1 - cos6U)(l - cos9 cd ) 



1 + 



(1 — cos 9 ac )(l — cos 8 
(1 - cos6> ac )(l - cosflfcrf) 
(1 - cos6> ad )(l - cos6> bc ) 

+ i27T 2 a3)5 i2 \n ac )5 i2 \n bd ) 



bd) 



x In 



+ 2 In 



(1 - cos0 ad )(l - cos6> fec ) 
(1 - cos6U)(l - cos8 cd ) 
(1 - cos8 ad )(l - cos0 6c ) 



(A.2) 



The eigenfunctions and the eigenvalues can be exactly mapped as well. 



B. The correspondence between boost and dilatation 

The correspondence of the gluon distribution has an interesting property under boost. One 
can see from ( |2.5D that the dilatation in the x-coordinates, 

x » _> x '» = \ x » (B.l) 
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dilatation 




boost 



Figure 5: Boosting the qq pair in e + e annihilation corresponds to changing the size of the dipole 
in high energy scattering. 

translates into the boost in the y 3 direction 

y+ = \ y '+, y- = ^, y = y'- (b.2) 

This transformation relates the gluon distribution in the boosted frame to that generated by 
a squeezed dipole. For simplicity, assume that the initial quark-antiquark pair is oriented 
to the y 1 -axis (see fig. ||). A is related to the velocity v of the boost via 

A = — - — - , (B.3) 
7(1 + «) 

where 7 = l/Vl — v 2 as usual. Since v = cos 6 q in the new frame where 9 q is the angle of 
the quark jet, one has 

1 sin#„ , j . 

A = r = q -7T = « , (B.4) 

7(1 +v) l + cosfl, 1 h V 1 

consistently with ( |B.1| ). 

Under the dilatation (|B.1|), the energy distribution (|2.6D transforms as 



S'(x') = ±E{x) (B.5) 



The distribution in the boosted frame is 

\3 



Using the relation 



(l + costrj^ X d (l + cos 0' y 



v + cos _. 

cos6' = - (B.7) 

1 + v cos 9 y ' 
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one finds 



1 



7 3 (1 — v cos #') 3 

Similarly the number (or the charge) distribution transforms as 

1 dN[n(Q')} 



dN 



7 2 (1 — v cos 9') 2 



dn 



(B. 



(B.9) 
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